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ABSTRACT 


This  article  presents  the  vibration  of  a single-layered  graphene  sheet  resting  on  an  elastic  foundation  by  using 
nonlocal  first-order  shear  deformation  theory.  Equations  of  motion  for  a simply-supported  graphene  sheet  are  obtained  via 
a nonlocal  shear  deformation  theory.  Effects  of  nonlocal  parameter  as  well  as  length  of  graphene  sheet,  mode  numbers, 
three-parameter  of  foundation  and  thermal  parameter  are  discussed.  A comparison  example  is  presented  to  show  the 


accuracy  of  the  present  results. 
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INTRODUCTION 

Vibration  analysis  of  nanoplates  using  the  nonlocal  theory  of  elasticity  [1-4]  is  the  main  subject  of  research  works 
in  recent  years.  The  single-layered  graphene  sheet  (SLGS),  or  double -layered  graphene  sheet  (DLGS),  or  multi-layered 
grapheme  sheet  (MLGS)  has  high  resistance  and  unique  properties.  So,  all  of  them  are  used  for  manufacturing  of  many 
devices  such  as  oscillators,  clocks  and  sensor  devices.  Sakhaee-Pour  et  al.  [5]  presented  vibration  analysis  of  SLGSs  by 
using  a molecular  structural  mechanics  method.  Pradhan  and  Phadikar  [6]  studied  the  nonlocal  vibration  of  SL  and  DL 
nano-plates  via  classical  plate  theory  (CPT)  and  first-order  shear  deformation  plate  theory  (FSPT).  Ansari  et  al.  [7] 
investigated  the  vibrational  behavior  of  a SLGS  via  FSPT  and  differential  equations  are  solved  by  using  generalized 
differential  quadrature  method  for  various  boundary  conditions.  Pradhan  and  Kumar  [8]  presented  the  vibration  analyses  of 
an  orthotropic  SLGS  using  the  CPT  and  solved  the  governing  equations  of  motion  by  using  differential  quadrature  method. 
Satish  et  al.  [9]  presented  thermal  vibration  analyses  of  orthotropic  nanoplates  based  on  nonlocal  continuum  mechanics  for 
small  scale  effects.  Shen  et  al.  [10]  used  Galerkin’s  method  to  present  the  vibration  analysis  of  a SLGS  based 
nano-mechanical  sensor  via  nonlocal  Kirchhoff  plate  theory. 

Most  of  nanostructures  are  resting  on  two-parameter  elastic  foundation.  Ansari  et  al.  [11,  12]  studied  the  vibration 
analysis  of  a MLGS  using  the  FSPT  according  to  Winkler-type  foundation.  Murmu  and  Adhikari  [13]  investigated  the 
nonlocal  vibration  of  bonded  double  nanoplate  systems  according  to  Winkler-type  foundation.  Wang  et  al.  [14]  used  the 
nonlocal  theory  to  derive  the  nonlinear  governing  equations  for  double -layered  nanoplates  subjected  to  four  different 
boundary  conditions  according  to  Winkler-type  foundation.  Behfar  and  Naghdabadi  [15]  presented  nano  scale  vibration 
analysis  of  a MLGS  embedded  in  elastic  medium.  Chien  et  al.  [16]  investigated  nonlinear  vibration  of  laminated  plates 
resting  on  a nonlinear  elastic  medium.  Liew  et  al.  [17]  proposed  a continuum-based  plate  model  to  study  the  vibration 
behavior  of  MLGSs  that  are  embedded  in  an  elastic  matrix. 
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Pradhan  and  Murmu  [18]  employed  nonlocal  plate  theory  and  used  DQM  for  vibration  analyses  of  nano-SLGSs 
embedded  in  elastic  medium.  Pradhan  and  Kumar  [19]  presented  vibration  analysis  of  orthotropic  SLGS  embedded  in 
Pasternak’s  elastic  medium. 


The  natural  vibration  frequency  of  a SLGS  resting  on  a two-parameter  Pasternak’s  foundation  [20-24] 
is  investigated  in  the  present  article.  The  nonlocal  elasticity  theory  via  the  first-order  shear  deformation  theory  is  presented. 
The  differential  governing  equations  are  derived  and  their  solution  is  analytically  presented  for  a simply-supported  SLGS. 
The  effects  of  nonlocal  index,  two-parameter  foundation  on  the  natural  vibration  frequencies  are  illustrated.  A comparison 
with  the  literature  is  presented  and  benchmark  results  are  plotted  for  sensing  the  effect  of  all  used  parameters. 

BASIC  EQUATIONS 


Let  us  consider  a SLGS  resting  on  elastic  foundation  and  subjected  to  uniform  load.  The  SLGS  is  of  length  a, 
width  b and  uniform  thickness  h as  shown  in  Figure  1 . The  SLGS  is  made  of  a homogeneous  isotropic  and  linearly  elastic 
material  with  Young’s  modulus  of  elasticity  E,  Poisson’s  ratio  v,  shear  modulus  G and  material  density  p. 


Figure  1:  A Continuum  Plate  Model  of  a Single-Layered 
Graphene  Sheet  Resting  on  Elastic  Foundation 


The  displacements  u;  of  a SLGS  according  to  the  FSPT  are  expressed  as 
ux {x, y, z,t)  —u ( x , y,t)  + zip ( x , y,  t), 

u2(x,y,z,  t)  = v(x,y,  t)  + z (p(x,y,  t),  (1) 

u3  ( x , y,  z,t)  = w ( x , y,  t ), 


in  which  u1,  u2,  and  u3  are  the  displacements  in  the  x,  y,  and  z directions;  u and  v are  in-plane  displacements 
while  w denotes  transverse  displacement  (deflection);  ip  and  0 are  rotational  displacement  about  y-  and  x-axes, 
respectively. 


The  strain-displacement  equations  of  elasticity  are  given  by 
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The  constitutive  equations  of  an  isotropic  SLGS  according  to  the  nonlocal  elasticity  can  be  expressed  as 
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where  <f  = (ae0)2  represents  the  nonlocal  index  in  which  the  parameter  a denotes  an  internal  characteristic  length 
and  e0  denotes  a material  constant  determined  experimentally. 

Hamilton’s  principle  is  used  here  to  obtain  the  equations  of  motion  as 

/t  [ fjj  (p  gt 2 + CTxxSexx  + (TyySEyy  + ^ Xy  ^ Yxy  ®xz^Vxz  4 Gyz^Yyz}  ^12  dz 


— jn  (q  — Rf^Sw  dnj  dt  = 0, 

where  q is  the  transverse  distributed  load  and  Rf  is  thetwo-parameter  Pasternak’s  foundation 
Rf  = (/^  -K2V2)w, 


(4) 


(5) 


in  which  V2  is  the  Laplacian,  K1  is  the  linear  Winkler’s  modulus  and  K2  denotes  Pasternak’s  (shear)  foundation 


modulus. 


Integrating  the  displacement  gradients  in  £j;-  by  parts  in  Eq.  (4)  to  obtain  the  governing  equations  in  the  form 
(setting  q — 0) 
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dx  dy 

where  Ij  denote  mass  moments  of  inertia  defined  by 

h = S-h/z  PzJ'Az- 


(7) 


and  Ntj,  Mtj , and  QL  denote  basic  components  of  stress  resultants,  stress  couples  and  shear  stress  resultants. 
They  can  be  obtained  by  integrating  Eq.  (3)  over  the  thickness  of  the  SLGS  as 


{%  Mip  <?;}  = f_h2/2{(Jij,  zatj,  kaiz]dz , i,j  — x,  y, 


(8) 


k is  the  transverse  shear  correction  factor.  After  using  Eqs.  the  above  stress  resultants  can  be  written  in  terms  of 
the  displacements  as 
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where  D = 


— — represents  bending  rigidity  of  the  SLGS. 


The  substitution  of  Eqs.  (9)-(ll)  into  Eqs.  (6)  yields  the  following  nonlocal  partial  differential  equations  of  motion 
in  terms  of  displacements 
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NONLOCAL  VIBRATION  FREQUENCIES 

The  determination  of  natural  frequencies  is  of  fundamental  importance  in  design  of  many  nano-structures. 
The  assumed  form  of  displacement  components  is  expressed  as 

({u,  xp}\  f {hu*,  xp*}  cos(Xmx)  sin(/rny)'J 

| w | = •!  hw*  sin(/l)nx)  sin(//ny)  (17) 

l{v,  cp})  [{hv*,  cp*}  sin(Amx)  cos (/rny) J 

where  a>  represents  the  natural  frequency  for  the  SLGS;  u* , v*,  w* , xp*  and  cp*  are  arbitrary  parameters; 

Xm  = mn/a  and  fin  = nn/b  in  which  m and  n are  the  wave  numbers;  and  i = V— 1. 

The  substitution  of  the  solution  given  above  into  equations  of  motion,  for  constant  density,  yields: 


D (x2m  + i^/12)  u*  + ^DXmfinv*  = ^<o2[l  + tal  + H2n)]u *,  (18) 

^DXmiinu*  + D(^X2m  + fi2)  v*  = ^tu2[l  + fal  + a£)]v*.  (19) 

kGh(X2n  + ^ )w*  + [1  + f(A2n  + /r?2„)][(fi  + (A2n  + 

+kGXmxp*  + kGfincp * = phoj2[l  + + //m)]w*<  (20) 

kGh2 Xmw*  + [D  (A2n  + -j-Pn)  + feGh]  xp*  + ^ DXmfinCp * = p-£ tn2[l  + fU2n  + p2)]0‘,  (21) 

kGh2finW*  + ^DXmHn r + [d  + Mn)  + fcG/l]  0*  = ^ tU2(l  + fU2„  + p2)]0*.  (22) 

The  above  equations  of  motion  may  be  written  in  terms  of  displacement  parameters  w*,  xp *,  and  <p*  only  as 
([P]-tu2[R])W  = {0],  (23) 
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where  {A'}  = {w* ,ip* ,(p*}T  is  the  solution  vector.  The  non-vanishing  elements  of  the  symmetric  matrices  | P|  and 
[R]  are  expressed  as: 

Pn  = kGh3(Al i + p2)  + h2[  1 + <f(Am  + /tm)][^l  + (^m  + 

P12  = kGh2Amxp*, 

P22  = D (a2^  + — Mn)  + kGh, 

P13  = kGh2pncp*, 

^23  2 P^mPn> 

P33  = D + n2^  + kGh, 

RX1  = ph3w2[l + +/4)], 

P22  =P33  = P-§-a>2[l+^A2m+p2)].  (24) 

The  frequency  equation  for  the  SLGS  is  given  by  setting  | \P\  — a>2 [R]  | = 0 to  get 

A3(o6  — A2co4  + A-lCO2  — A0  — 0,  (25) 

in  which 

^0  = ^22(^11^33  — Pi  3)  ~ P\2  (^12^33  — P13P23 ) — P23  (P\  1 P23  ~ P12P13)' 

A-l  — Rll  [P22P33  ~ P23)  + ^22^1^33  — P\3 ) + P22  (P11P22  ~ P12)’ 

A2  ~ [^11^22  + (P22  + ^,33)^ll]^22< 

•4  3 = P 11^22-  (26) 

NUMERICAL  RESULTS  AND  DISCUSSIONS 

In  this  section  the  properties  of  graphene  sheet  are  considered  as  [6]:  modulus  of  elasticity  E — 1.02  TPa, 
Poisson’s  ratio  v = 0.16  and  material  density  p — 2250  kg/m3.  It  is  assumed,  except  otherwise  stated, 
that  a — b — 10  nm  and  h — 0.34  nm.  Here,  the  fundamental  vibration  frequency  (m  = n = 1)  is  firstly  compared  with 
the  corresponding  ones  in  the  literature  [6,  8,  10,  19].  For  this  purpose,  the  frequency  ratio  (a>  — o)NL/o)L)  is  considered 
where  o> N L represents  the  ftrequency  calculated  using  nonlocal  theory  while  o>L  represents  the  Frequency  calculated  using 
local  theory.  The  comparison  of  vibration  frequency  ratio  of  nonlocal  square  SLGS  is  reported  in  Table  1.  An  excelent 
agreement  is  appeared  and  the  present  ratio  is  the  same  as  those  presented  in  the  literature. 

Now,  let  us  consider  additional  examples  to  put  into  evidence  the  effect  of  length  a,  the  foundation  parameters 
/fj  and  K2  and  the  nonlocal  index  on  the  vibration  of  the  present  SLGS.  It  is  to  be  noted  that,  we  can  get  the  local  thermal 
vibration  of  the  SLGSs  by  setting  = 0 in  the  preceeding  equations. 
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Table  1:  Comparison  of  Fundamental  Vibration  Frequency  Ratio  co  of 
Nonlocal  Square  SLGS  (a  — b — 10  nm  and  h — 0.  34  nm) 


? (nm2) 

Present 

[6] 

[8] 

[10] 

[19] 

0 

1 

1 

1 

1 

1 

1 

0.91386 

0.9139 

0.9139 

0.9139 

0.9139 

2 

0.84673 

0.8467 

0.8466 

0.8467 

0.8468 

3 

0.79251 

0.7925 

0.7926 

0.7925 

0.7926 

Figure  2:  Fundamental  Frequency  Ratio  of  a SLGS  Vs  Nonlocal 
Index  and  its  Length  ( b — lOnm,  h — 0. 34nm,  K1  = 0,  K2  — 0) 


Figure  3:  Natural  Frequency  Ratio  of  a Square  SLGS  Vs  Nonlocal 
Index  (a  = b — lONm,  h — 0. 34Nm,  m = 1,  K±  = 0,  K2  — 0) 


Figure  4:  Fundamental  Frequency  Ratio  of  a Square  SLGS  Vs  Nonlocal 
Index  and  Winkler’s  Parameter  (h  — 0. 34nm,  a = b = lOnm,  K2  — 5) 
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Figure  5:  Fundamental  Frequency  Ratio  of  a Square  SLGS  Vs  Nonlocal 
Index  and  Pasternak’s  Parameter  ( h — 0. 34  nm,  a = b — 10  nm,  K1  — 10) 


Figure  2 shows  the  fundamental  frequency  ratio  a>  vs  nonlocal  index  <f  for  different  values  of  the  length  a of  the 
SLGS  without  elastic  foundation  (Kt  = K2  — 0).  The  frequency  ratio  63  increases  as  a increases  and  as  <f  decreases. 
The  maximum  fundamental  frequency  ratio  occurs  when  f = 0.  Figure  3 shows  the  natural  (m  = 1 and  n = 1, 2,3,4) 
frequency  ratios  63  vs  nonlocal  index  f of  the  SLGS  without  elastic  foundation  (K±  — K2  — 0).  The  natural  frequency  ratio 
63  increases  as  n decreases. 


Figures  4 and  5 show  the  fundamental  frequency  ratio  S of  a square  SLGS  vs  nonlocal  index  for  different 
values  of  Winkler’s  parameter  K±  and  Pasternak’s  parameter  K2,  respectively.  Once  again,  the  fundamental  frequency  ratio 
ay  increases  as  ( decreases.  The  maximum  frequency  ratio  occurs  when  = 0 and  for  all  values  of  K1  and  K2.  The 
frequency  ratio  is  increasing  with  the  increase  of  Winkler’s  parameter  K1  and  Pasternak’s  parameter  K2. 


Figure  6:  Fundamental  Frequency  Ratio  of  a Square  SLGS  Vs 
Pasternak’s  Parameter  ( h = 0.  34  nm,  a — b = 10  nm,  £ =1) 


Figure  7:  Fundamental  Frequency  Ratio  of  a Square  SLGS  Vs 
Winkler’s  Parameter  ( h — 0. 34  Nm,  a — b — 10  Nm,  ^ =1) 
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The  plots  of  fundamental  frequency  ratio  ai  of  a square  SLGS  vs  the  elastic  foundation  parameters  Kt  and  K2  are 
appeared  in  Figures  6 and  7.  The  frequency  ratio  is  independent  of  K2  for  small  values  of  K1 . For  constant  values  of  K1  the 
frequency  ratio  is  increasing  as  K2  increasing  (Figure  6).  Also,  for  constant  values  of  K2  the  frequency  ratio  is  rapidly 
increasing  as  K1  increasing  (Figure  7),  especially  in  the  range  8 < K,  < 14. 

CONCLUSIONS 

The  nonlocal  vibration  frequencies  of  a single-layered  graphene  sheet  resting  on  two-parameter  Pasternak’s 
foundation  are  investigated.  The  local-to-nonlocal  vibration  frequency  ratios  are  compared  well  with  the  corresponding 
results  in  the  literature.  The  fundamental  and  natural  frequency  ratios  are  presented  to  serve  as  benchmarks  for  future 
comparisons.  It  can  be  observed  some  novel  phenomena  from  the  discussion  of  the  results.  The  results  are  very  sensitive  to 
the  variation  of  nonlocal  index.  The  inclusion  of  two-parameter  Pasternak’s  foundation  is  very  pronounced.  The  maximum 
frequency  ratio  occurs  when  neglecting  the  nonlocal  index. 
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